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The Hubble constant from galaxy lenses: impacts of triaxiality and 
model degeneracies 
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ABSTRACT 

The Hubble constant can be constrained using the time delays between multiple images 
of gravitationally lensed sources. In some notable cases, typical lensing analyses assuming 
isothermal galaxy density profiles produce low values for the Hubble constant, inconsistent 
with the result of the HST Key Project (72 ± 8 kms -1 Mpc -1 ). Possible systematics in the 
values of the Hubble constant derived from galaxy lensing systems can result from a number 
of factors, e.g. neglect of environmental effects, assumption of isothermality, or contamination 
by line-of-sight structures. One additional potentially important factor is the triaxial structure 
of the lensing galaxy halo; most lens models account for halo shape simply by perturbing the 
projected spherical lensing potential, an approximation that is often necessary but that is in- 
adequate at the levels of triaxiality predicted in the CDM paradigm. To quantify the potential 
error introduced by this assumption in estimates of the Hubble parameter, we strongly lens 
a distant galaxy through a sample of triaxial softened isothermal halos and use an MCMC 
method to constrain the lensing halo profile and the Hubble parameter from the resulting mul- 
tiple image systems. We explore the major degeneracies between the Hubble parameter and 
several parameters of the lensing model, finding that without a way to accurately break these 
degeneracies accurate estimates of the Hubble parameter are not possible. Crucially, we find 
that triaxiality does not significantly bias estimates of the Hubble constant, and offer an ana- 
lytic explanation for this behaviour in the case of isothermal profiles. Neglected triaxial halo 
shape cannot contribute to the low Hubble constant values derived in a number of galaxy lens 
systems. 

Key words: Gravitational lensing - Galaxies: fundamental parameters - galaxies: halos - 
cosmology: cosmological parameters - cosmology: theory - dark matter 



1 INTRODUCTION 

Gravitationally lensed, multiply imaged radio source time de- 
lays are used to constrain the Hubble constant. Lensing analy- 
ses produce a range of values, some of which are significantly 
lower th an the r esult of the HST Key Project (72 ± 8 kms" 
Mpc -1 ); lOguril J2007l) offers an excellent summary of the cur- 
rent state of the field. Further, although there is a general agree- 
ment upon is othermality from a number of differently motivated 
studies (e.g. IrTx et al. I 1 19971: iRomanowskv & Kochanekl 1 19991: 



Gerhard et all 120011: i Treu & Koopm ans] 1 20021 : iKoopmans et al.l 
2006l : ISaha et alJl2006l ; lDobke & Kindl2006h . modelling of a num- 
ber of systems (e.g. PG1115+080, B1600+434, HE2149-2745, 
and SBS 1520+530) has shown that these pref er more centrally 
concentrated, non-isotherma l, density profiles (Kochanek 120021 : 
iKochanek & Sche chter 2004). The slope of the lensing density pro- 
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file and the value of Hubble constant are heavily degenerate, and 
the number of constraints from galaxy lens systems small, making 
it unclear whether the inconsistencies in slope and Hubble constant 
values indicate flaws in the measurement of both or only one of the 
two quantities. 

While some recent statistical studies (e.g. ISaha et all d2006h ) 
find values consistent with the consensus value, some unexplained 
inconsistencies remain among the varied results across the field, 
and whether these stem primarily from external environmental ef- 
fects, intrinsic physical attributes of the lensing galaxy or universe, 
or differences in modeling techniques is still unclear. A number of 
factors are known to influen ce the derivation of the Hub ble constant 
from galaxy lens systems. iKeeton & Zabludoff (2004) have shown 
that lens models of four-image systems overestimate the Hubble 
constant by up to 15% when neglecting the effect of the contribu- 
tion from other group members on the potential of the lens galaxy, 
with an even higher discrepancy for two-image systems. Addition- 
ally, tidal stripping and shocking has been shown to induce tran- 
sient fluctuations in the inner density profile of lens galaxies present 
in groups, and the resultant profile steepening can also contribute 
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to erroneous Hubble constant v alues if isothermality is assumed 
( Dob ke. King & Fellhaueif 2007). Outside the group environment, 
line-of si ght structures can also provide contributions to the lens 
potential dMomcheva e t al. 2006). These effects can result in either 
an over or underestimate of the Hubble constant by virtue of a con- 
tributed increase or decrease in the density profile slope. 

The intrinsic shape of the lensing galaxy has also been sus- 
pected to ha ve a direct bearing on Hubble constant derivation. 
lOguril d2007l) recently demonstrated the importance of substruc- 
tures and external perturbations in galaxy lensin g; furthermore, in 
an inv estigation with models of 35 galaxy lenses ISaha & Williams! 
(2006) proposed that shape-modelling degeneracies (e.g. caused by 
triaxiality) could also contribute to changes in the time delays and 
hence the derived Hubble constant values. Triaxiality is of par- 
ticular interest because the dark matter halos of ACDM are pre- 
dicted to be fully t riaxial with axis ratios in the mass distribution 
as low as 0.4 (e.g. iBett et alJkoOTt iMaccio et ai1l2q07l). a predic- 
tion w hich observations have loosely confirmed (e.g. lBak & StatleJ 
2000). Moreover, it has been shown that the triaxiality of CDM ha- 
los can effect the overall lensing probabilities and relative number 
of diffe rent image configurations (double, quadr uple, naked cusp 
lenses) dOguri & Keeto n 2004; Rozo et aI1l2007l) . Triaxial studies 
in galaxy cluster weak lensing have further revealed that neglect- 
ing halo shape can result in sign ificant over an d underestimates 
of cluster concentration and mass dCorless & King||2007l) . It seems 
apparent then that triaxiality can affect a number of critical observ- 
ables and derived quantities in both strong and weak lensing and 
as such cannot be ignored. However, due to the small number of 
constraints available in galaxy lensing analyses, halo shape is typ- 
ically addressed only via an elliptical perturbation to the projected 
spherical lensing potential, ignoring altogether halo structure along 
the line-of-sight. 

In this paper we investigate what errors this often necessary 
neglect of triaxiality in galaxy lensing analyses introduces in best 
fit Hubble parameter values. We fit an elliptical isothermal mass 
model to lensing systems of highly (as predicted in simulations) 
triaxial isothermal galaxy lenses to investigate the maximum dis- 
crepancies that can arise from fitting a simplified model to data that 
in fact originates from a fully triaxial lens. We ask: to what ex- 
tent can galaxy triaxiality explain the discrepancies between some 
lensing-derived values for the Hubble constant and those of other 
methods? The outline of the paper is as follows: in §2 we intro- 
duce the softened triaxial isothermal model and its lensing prop- 
erties, §3 goes on to present our analysis method and key results, 
and §4 draws some conclusions based upon our findings. A stan- 
dard ACDM cosmology with H = 72 km/s/Mpc, Q, m = 0.3, 
Q.a = 0.7 is employed throughout. 



2 LENSING BY A SOFTENED TRIAXIAL ISOTHERMAL 
HALO 

To generate a full parameterisation for a soften ed triaxial isother - 
mal halo we follow the procedure giv en first inljing & Sutol ( 2002) 
and implemented for weak lensing in lCorless & K ing ( 2007}) for a 
triaxial NFW. We first generalise the spherical softened isothermal 
profile to obtain a density profile 
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where 9e is an effective Einstein radius, S a triaxial core radius, 
a triaxial radius 
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a/c and b/c the minor:major and intermediate:major axis ratios, 
respectively, and Di a the angular diameter distance between the 
lens and source and D s that between the observer and source. 
The mass contained within radius R is 
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The virial mass is defined as that mass contained within an ellipsoid 
containing a mean density 200 times the critical density p c (z) at the 
redshift of the halo, M200 = ^^-abR^ooPciz); combining that 
definition with ([3} gives an expression for 9e as a function of the 
virial radius R200 '■ 
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2.1 Lensing Properties 

The full derivation of the le n sing properties of a triaxial halo is 
given by lOguri. Lee, & Sutol d2003l) . and we summarise some of 
that work here. The triaxial halo is projected onto the plane of the 
sky to find its projected elliptical isodensity contours as a function 
of the halo's axis ratios and orientation angles (9, (j>) with respect 
to the the observer's line-of-sight. The axis ratio q of the elliptical 
density contours is given by 
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the elliptical radius by 
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and the orientation angle Q, on the sky by 
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General expressions for the l ensing potentia l of a softened el- 
liptical isothermal halo are given in lKeetonl JioOl): we modify them 
to incorporate the geometry of the halo projection to establish a di- 
rect connection between a 3D halo and its lensing potential. The 
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lensing potential is 



E qS\n J (r + S)* + (1 - q 2 ) X* 



where 

r 2 = q 2 (s 2 + x' 2 )+y 2 , 



X 



@Eq , -i 
: tan 



V^v 2 



VT^q 2 * 

r + s 



T + q 2 S 



, and 



qx qx 



(14) 

(15) 
(16) 

(17) 
(18) 



where X, Y, and S are all angular distances measured in the im- 
age plane. These analytic formulations were confirmed using the 
independent numerical methods for calculating lensing properties 
of elliptical density profiles by integ rating over those of their spher- 
ical counterp arts as given inlKeetonl fcOOll) and applied to a triaxial 
NFW halo in lCorless & Kind j2007l) . as well as by brute force pro- 
jection of the triaxial halo. 



2.2 Time delays and the Hubble parameter 

The time delays and image positions of strongly lensed systems 
are completely constrained by th e light-travel time equation (e.g. 
ISchneider. Ehlers, & Falcol ( 1 19921) ) 
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respectively, where 9 and j3 are angular positions in the lens and 
source planes respectively, ip is the lensing potential, and z\ is the 
redshift of the lens. Examining these with some care shows that it 
is the image positions, and not the time delays between them, that 
set the normalisation of the model potential: when matching im- 
age positions one matches the stationary points of the time-delay 
surface, thus mandating the model satisfy only the lens equation 
9 — /3 = SJip. This is observationally independent of the Hubble 
parameter, for while these are all angular values, with 9 = X / Di , 
P = X /D a , and the physical distances X are not known and can- 
not therefore be used as a constraint on D s and Di . Once the shape 
and normalisation of the potential is set by the image position con- 
straints via the first derivative of the time delay surface, the nor- 
malisation factor DiD s /Di s is crucially important in setting the 
time delays between the multiple images, and it is here the Hubble 
parameter can be constrained. 
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Figure 1. Mean marginalized parameters for oblate and prolate halos with 
(a = 0.6, b = 1.0) and (a = 0.6, b = 0.6) for fixed slope parameters 
r\ = {0.9, 1.0, 1.1} rotated through 90° from an orientation of projected 
circular symmetry at 0° to that of maximum visible ellipticity. The dotted 
line shows the true parameter value of the underlying lens. 



3 RESULTS & DISCUSSION 
3.1 Simulations & Fitting Method 

To generically model the impact of unaccounted for triaxial struc- 
ture in isothermal galaxies we use a Markov Chain Monte Carlo 
(MCMC) sampling method to explore the posterior probability dis- 
tribution obtained when fitting an elliptical power law model to a 
triaxial isothermal lens. First, simulated multiple-image lens sys- 
tems are generated by lensing through a triaxial isothermal halo 
according to Eq. [14] and the lens equation. These simulated ob- 
servations of image positions and time delays are then fit using 
a general 2D elliptical power law model governed by seven free 
parameters. Such models cannot be constrained using two-image 
multiple-image systems, as these provide only 5 constraints (2 im- 
age positions vectors + 1 time delay); 4-image systems fare better, 
providing 1 1 constraints (4 image position vectors + 3 time delays). 
However, even in 4-image galaxy lens systems, the data is typically 
able to meaningfully constrain only one of either the slope of the 
underlying galaxy potential or the Hubble parameter, due to an oft- 
noted strong degeneracy between these two key parameters. 

Several papers have explored constraints on H from multi- 
ple image time delays, and degeneracies with parameters describ- 
ing th e lens model, in particu lar t he radial profile. Early w ork in- 
cludeslGorenstein et aildl988l) and lRefsdal & Surdeil dl994l) . Later, 
IWitt et alJd2000h consider analytic time delays and estimates of H a 
from a general family of lensing potentials which allows for an- 
gular structure. Their treatment w as extended and generalised by 
IWucknitzl ( l2002l) . lKochanek1 d2002l) proposes a semianalytic model 
to understand observed time delays, isolating the important phys- 
ical properti es of the lens models. For example, building on the 
treatment of iGorenstein et all dl988h . he obtains scaling solutions 
for H for general lenses which depend on the mean convergence 
inside annuli bounded by the multiple images. Due to the strong de- 
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Figure 2. The posterior probability distribution of a power law model with fixed isothermal slope r] = 1 fit to an isothermal triaxial prolate halo with axis 
ratios (a = 0.6, b = 0.6) at 80°, giving elliptical isodensity contours of axis ratio 0.61. The parameter giving the orientation angle of the isopotential contours 
9 is omitted for conciseness, as it shows no degeneracy with any other parameters and is always returned tightly constrained to the true value. The shading 
(dotted lines) gives the 2D (ID) likelihood distribution. 



generacy, in most of these and other analyses some assumption is 
often made regarding the shape of the galaxy density profile and the 
Hubble parameter then obtained under that assumption (though si- 
multane ous fits have been succesfully carried out in statistical stud- 
ies, e.g. ISaha et aT] J2006)). Most often the galaxy halo is taken 
to be isothermal with density profile slope p(r) oc r~ 2 , a well- 
motivated simplificati on derived from previo us studies of galaxy 
halo profiles (see eg. Koopmans et al.l J2006l) ). However, for any 
one galaxy, it is quite unlikely that the slope of the profile will be 
exactly isothermal. Given the strong degeneracy between the Hub- 
ble parameter and the density profile slope, such systematic errors 
may be quite important. 

To avoid having to make such a potentially error-inducing as- 
sumption, we adopt an MCMC fitting method that is capable of 
exploring the full posterior probability distribution and all its de- 
generacies. MCMC methods employ a "guided" random walk that 
returns a sample of points representative of the posterior probability 
distribution; the probability of a certain region of parameter space 
containing the true model is directly proportional to the density 



of points sampled in that region. From the distribution of sample 
points the full posterior probability distribution is obtained, which 
is easily and direcly marginalized over to obtain fully marginal- 
ized most-probable parameter estimates for all parameters. These 
most-probable parameter values reflect the full shape of the pos- 
terior probability distribution without assumption about the error 
distribution, and with full knowledge of the prior. Such methods 
have exploded in popularity recently, and there a re several excel- 
lent re f erences describi ng the method in detail, e.g. lLewis & Bridle 
( 2002), MacKav ( 2003). Simply put, the sampler functions by step- 
ping through parameter space by taking random steps governed by 
a transfer function, usually a simple n-D Gaussian, where n is the 
number of parameters of the fitted model. If the randomly-chosen 
next step is to a point of higher probability than the current posi- 
tion, the step is taken. If the next step is to a point of lower proba- 
bility, the step is taken with probability p(new) /p(current). Thus, 
the MCMC sampler spends most of its time in high probability re- 
gions, but can move "downhill" to regions of lower probability in 
order to explore the entire space and sample all regions of high 
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probability. Crucially, this method is able to return a true represen- 
tation of the full posterior probability distribution, regardless of the 
complexities of that probability distribution, e.g., tight degenera- 
cies, multiple islands of high probability, or a very flat distribution 
due to poor constraints from the data. 

Thus, we are able to fully understand the impacts of fitting 
2D models to systems with 3D triaxial structure, without induc- 
ing errors by making other limiting assumptions. To implement the 
method, we must define the posterior probability function, defined 
in Bayesian statistics as 

p(%)p(t) 



p(n\6) 



p{6) 



(21) 



where p(9\n) is the likelihood of the data given the model param- 
eters (the standard likelihood), p(ir) is the prior probability dis- 
tribution for the model parameters (e.g. an outside constraint on 
the slope of the density profile), and p(9) is a normalising factor 
called the evidence, of great value in comparing models of differ- 
ent classes and parameter types, but expensive to calculate and un- 
necessary for the accurate exploration of the posterior distribution. 
The assignment of priors is often a controversial exercise, but in this 
case it is straightforward to assign weak, physically sensible priors 
that do not significantly affect our results. We define the likelihood 
function as a standard exponential function of the differences be- 
tween the model image positions and time delays and those of the 
"observed" data: 



C = n i=1 exp 
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n 



=1 exp 
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(22) 



and take the measurement errors in the image positions a x and time 
delays at to respectively be 0.7% and 14%, reasonable values for 
strong lensing analyses. 

In our MCMC sampler we employ a 7D two-sided Gaussian 
transfer function, and use the covariance matrix of an early run to 
sample in an optimised basis aligned with the major degeneracies 
of the posterior (of which there are several!). We tune the step sizes 
of the sampler to achieve an average acceptance rate of 1/3 in each 
basis direction, run three independent MCMC chains, started at 
randomly chosen positions in parameter space, for each lens sys- 
tem, and sample the distribution space until the standard var(chain 
mean)/mean(chain var) indicator is less than 0.2, indicating chain 
convergenc e. We utilise the Get Dist package from the standard 
CosmoMC dLewis &~B ridle 2002) distribution to calculate conver- 
gence statistics, parameter contours, and marginalized parameter 
estimates. The highly probable region of parameter space is of very 
small volume compared to the that of the prior parameter space, 
and we find it is a significant time-saver to first follow an MCMC 
chain with very large step sizes to find regions of high probabil- 
ity, then restart a chain at this randomly found, but advantageously 
chosen, starting point with a step size tuned to the size and shape 
of the posterior distribution. This does not upset the statistics of the 
MCMC method; it is equivalent to throwing away points from the 
very long burn-in phase that would be required if the step size were 
kept tuned to the small steps best for exploring the high probability 
region. 

In every case we fit a power-law lensing density profi le model 
to the simulated lensed images and time delays, given by iBarkanal 
( 1998) (here written with a factor of 2 modification to match a stan- 
dard elliptical isothermal model when r\ = 1), with lensing poten- 
tial 
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Figure 3. Mean marginalized parameters for oblate and prolate halos with 
a range of triaxial minor axis ratios a = {0.5, 0.6, 0.7, 0.8} at a rotation 
angle of 80° for fixed slope parameters r\ = {0.9, 1.0, 1.1}. The dotted 
lines show the true parameter values of the underlying lens. 
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where q p is the axis ratio of the isopotential ellipses. Note that 
here the ellipticity is in the 2D potential, as opposed to our triaxial 
model in which it is introduced in the 3D isodensity surfaces. At 
high ellipticities (e p > 0.3) where the projected axis ratio is small 
(q p < 0.54), near-isothermal elliptical lensing potential models 
break down, giving unphysical "dumbbe lF'-shaped projected mass 
distributions dKassiola & Kovnen 1 119931) ). However, in the regime 
of ellipticities we study here (q p > 0.75), this problem is insignifi- 
cant, and we choose this elliptical potential model as representative 
of those most often fit in lensing analyses and for its analytical sim- 
plicity. The potential is always less elliptical than the underlying 
mass distribution, such that the axis ratio q p is always greater than 
the axis ratio q of the isodensity contours. The underlying density 
profile steepens with decreasing r\, such that r\ — corresponds 
to a Hubble density profile that goes as r~ A . Thus, the seven free 
parameters are the profile slope parameter 77, the effective Einstein 
radius fee, the axis ratio of the model isopotential contours q p , the 
orientation angle of those contours 6, the source position compo- 
nents /3i and fa, and the Hubble parameter h. 

We begin our analysis with a simple case, similar in method to 
many typical lensing analyses, in which we set the slope parameter 
to a range of fixed values and carry out parameter fits, and exam- 
ine the impacts of triaxiality in that simplified framework. We then 
offer an analytic explanation for our numerical results. Finally, we 
examine the problem using the fully generalised power-law model 
in which the slope is allowed to vary freely, to complete our under- 
standing of the role of triaxiality and to better understand the larger 
degeneracies that govern the problem of Hubble constant determi- 
nation from strong lensing systems. 
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Figure 4. The posterior probability distribution of a general power law model fit to an isothermal triaxial prolate halo with axis ratios (a = 0.6, b = 0.6) 
oriented at 80°, giving elliptical isodensity contours of axis ratio 0.61. The parameter giving the orientation angle of the isopotential contours 9 is omitted for 
conciseness, as it shows no degeneracy with any other parameters and is always returned tightly constrained to the true value. The shading (dotted lines) gives 
the 2D (ID) likelihood distribution. 



3.2 A Fixed-slope model fit to oblate and prolate halos 

To begin to isolate the impacts of triaxiality in the problem of Hub- 
ble parameter determination from strong galaxy lensing systems, 
and to allow for direct comparisons to s tandard methods employed 
in many analyses (e.g. Koc hanek et al.l ( I2006al) ). we fix the density 
profile slope at one of three values and carry out an MCMC explo- 
ration of parameter space. We fix the slope of the underlying halo 
to either isothermal (77 = 1), steeper than isothermal (77 = 0.9), 
or shallower than isothermal (77 = 1.1), approximating the errors 



between the true and model slopes sure to sometimes occur if the 
slope is externally constrained to avoid the strong degeneracy be- 
tween it and the Hubble parameter. We then fit these models to 
lens systems of symmetric oblate and prolate triaxial lensing halos 
at z — 0.5 with a lpc core with axis ratios (a — 0.6,6 = 1.0) 
and (a = 0.6, b — 0.6), respectively, each containing a mass of 
1 x 10 11 Mq inside the triaxial effective Einstein radius Re = 
X 2 E /a 2 +Y E /b 2 +Z%=6 E D l . 
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Figure 5. Projections of the posterior probability distribution into the r\-h plane of a power law model fit to an isothermal triaxial prolate halo with axis ratios 
(a = 0.6, b = 0.6) oriented at 80°, giving elliptical isodensity contours of axis ratio 0.61 , shaded by the other important model parameters, respectively, bE, 
q p , 9, and /3i . Note the strong interdependence of the r\-h degeneracy on bp;, q p , and /3. 



3.2.1 Hubble parameter as a function of triaxial orientation 

We fit to each halo type at ten different orientations stepped from 
0° to 90°; at 90° each halo looks maximally elliptical with pro- 
jected isodensity axis ratio a while at 0° the projected profile is 
circular. When the projected ellipticity is small only 2 images are 
formed, and we discard these systems because they provide too few 
constraints to meaningfully study the Hubble parameter, leaving 
7 prolate systems and 5 oblate systems for study (there are more 
4-image prolate systems because prolate halos are more efficient 
lenses than oblate, such th at their threshold for 4- image systems 
is at a lower ellipticity: see ICorless & King] d2007h ). At each ori- 
entation the halos lens a source at z = 2 positioned inside the 
multiple-imaging region at /3 — {0.02", 0.03"}. The image posi- 
tions are found via the lens equation using a standard rootfinder, 
and the time delays between images calculated via the time de- 
lay surface (Eq. 1 19b; central images are discarded due to their high 
level of d emagnification, which makes them invisible in observa- 
tions (e.g. lWallington & Naravanl dl993h ). 

The results are plotted in Figure [TJ Crucially, when the slope 
of the fitted halo matches that of the underlying lens - when they 
are both isothermal - the Hubble parameter is accurately returned 
in every case! Thus, when the slopes are accurately matched in lens 
modelling, triaxiality induces no error in estimates of the Hubble 



parameter. The triaxial cases follow the familiar trends observed 

when there is a mismatch between the model slope and the true 

I 1 J — — k 

slope of the lens: as e.g. Kochanek (2002) has previously shown, 

fitting a model shallower than the true potential leads to underes- 
timates of the Hubble parameter, while overestimates are apparent 
when the fitted model has a steeper slope than the true lensing halo. 
The error in the returned Hubble parameter is worst when the visi- 
ble axis ratio is lowest and the lensing potential most elliptical. 



Figure[2]shows the posterior probability distribution when the 
power law model is fit to a prolate halo oriented at 80° with the 
slope set to match the true value of the isothermal lensing halo 
rj = 1. The orientation angle of the isopotential contours 9 is omit- 
ted for conciseness; it shows no degeneracies with any other param- 
eter and is very tightly constrained to the true value. Importantly, 
the only significant degeneracy visible is between the the compo- 
nents of the source position. Thus, when the slope is constrained 
accurately, a meaningful and accurate value for the Hubble param- 
eter can be derived from strong lens systems! However, the slope 
must indeed be accurately constrained; Figure[Tjshows clearly that 
a mismatch between the true slope value and the model slope value 
yields signficantly incorrect values for the Hubble parameter. 
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3.3 An analytic evaluation of the impact of triaxiality on Ho 

To better understand this rinding that Hubble parameter estima- 
tion is unaffected by the line-of-sight structure of triaxial halos, 
we look to the case of triaxial isothermal halos fitted with elliptical 
isothermal models, and more closely examine the expressions for 
the triaxial isothermal potential and its gradient, and their role in 
determing Ho in a strong lensing system. The potential of the soft- 
ened triaxial isothermal model was given in Eq |141 its derivatives 
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Figure 6. Marginalized mean best-fitting parameters for symmetric oblate 
and prolate halos of minor axis ratio a = 0.6 (filled symbols), as well as 
the maximum-likelihood parameters (unfilled), plotted as a function of the 
rotation angle of the lensing halo: at 0° the halos are circularly symmet- 
ric in projection, at 90° they show a maximum visible ellipticity of 0.6 in 
the projected isodensity contours. The dotted lines show the true parameter 
values of the underlying lensing halo. 



3.2.2 Hubble parameter as a function of triaxial axis ratio 

To complete this portion of the analysis, we seek to further clarify 
any impact of triaxiality by calculating the best-fitting parameter 
values for oblate and prolate halos with a range of triaxial minor 
axis ratios a = {0.5,0.6,0.7,0.8}, all oriented at 80° so as to 
have high projected ellipticities, at the same three fixed values of r\. 
The results are plotted in Figure [3] Again, when the density profile 
slope is mismatched to the true underlying value, increased pro- 
jected ellipticity exacerbates the error induced in the Hubble pa- 
rameter, with the best-fit value of h approaching slightly closer to 
the true value as the minor axis ratio a moves closer to 1 . That this 
trend occurs both for halos of the same fundamental shape oriented 
at different angles and for halo of different shape oriented at the 
same angle, and for both prolate and oblate halos, indicates it is 
purely associated with the visible ellipticity, with no dependence 
on the (very different) 3D structure of the various halos. 

We return to the point that there is essentially no difference in 
the value of the recovered Hubble parameter if the halo is oblate or 
prolate. As noted, the geometry of these two cases is quite different. 
When oriented at 80° a prolate halo is like a rugby ball being looked 
at edge on - its long axis is visible in the plane of the sky and 
there is relatively little mass along the line of sight. Conversely, the 
oblate halo is like a pancake looked at edge on, with its one short 
axis in the plain of the sky and a large amount of mass along the 
line of sight. This difference in projected mass is reflected in the 
potential normalisation &£. However, the Hubble parameter is in 
no way effected by this significant difference in 3D structure! The 
line-of-sight structure of the lensing halo appears unimportant to 
the Hubble parameter problem. 
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Examining Eqs. ll4|[24l and 1 25 1 we see that, when S is small, 
the triaxial structure of the halo manifests itself only through mul- 
tiplicative prefactors. This is not a generic behaviour for all poten- 
tials; it is true of this model because, when S is small, all terms 
carry factors of the same power of qx, the normalising ratio that 
carries information about the length of the axis hidden along the 
line of sight. All terms in the potential that are important when 
S is small are proportional to X or y, giving each term an over- 
all normalisation of qx, and similarly all important terms in its 
derivative are dimensionless in lengths giving each term a nor- 
malisation of qx- When S is large, or were there any terms that 
carried unequal powers of qx, this simple multiplicative renormal- 
isation due to triaxial structure would no longer hold. However, 
given that there is good eviden ce that galactic cores are very small 
dWallington & NaravarJ ( ll993l) ) and that galaxy density profiles are 
generally isothermal over the central region to which strong lensing 
is sensitive, this straightforward analysis indicates that triaxiality 
will not independently lead to any significant change in the shape 
of the lensing potential and thus creates no significant errors in the 
measurement of the Hubble parameter, as indicated in our simula- 
tions. The multiplicative factors introduced by the triaxial structure 
of the lensing halo are apparent in the varied values of the normal- 
isation of the potential 6e returned for different orientations of the 
triaxial halos. 



3.4 A general model fit to symmetric oblate and prolate halos 

To complete our understanding of the problem of Hubble parame- 
ter estimation from triaxial lens systems, we now look to the fully 
general case in which the elliptical power law model is fit with both 
Hubble and slope parameters left free to vary. 
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3.4.1 Hubble parameter as a function oftriaxial orientation: 
most-probable parameters 

We fit this general model to the same family of halos described in 
Section 13.2.1 1 The complete posterior probability distribution for 
one of the halos, the prolate halo oriented at 80° such that the ellip- 
ticity is approaching the maximum, is shown in Figure [4] again 
absent the well-constrained and independent orientation angle 9 
to make the plot clearer. The well-known degeneracy between the 
Hubble parameter h and the slope parameter rj is clearly visible, 
but interestingly, is not the only strong degeneracy present. Both 
the Hubble parameter and the density profile slope are also strongly 
degenerate with the normalisation of the potential 6b, and with the 
source position ft. 

The regions of highest likelihood (indicated by the shad- 
ing) are centred on the true parameter values; however, the most- 
probable parameter values obtained by marginalizing over all other 
parameters give consistently high values for the Hubble parame- 
ter and low values for the profile slope, indicating that while the 
true parameter values give a slightly higher likelihood value than 
others, there are other likely parameter value combinations includ- 
ing low values for r\ and high values for h that occupy signifi- 
cantly larger volumes in parameter space. Were the errors in this 
problem Gaussian, the maximum-likelihood parameter estimates 
should correspond to the most-probable marginalized parameter 
value; however, as seen in Figure [4] the posterior probability dis- 
tribution is highly non-Gaussian, leading to the separation of the 
point of maximum likelihood from the point of maximum proba- 
bility. The marginalized most-probable parameter values represent 
the true 'best' model; thus, left unbroken, the degeneracies of the 
problem produce real errors in any lensing result. 

To better picture the way these degeneracies function, Fig- 
ure [5] plots the posterior probability distribution in the r/-h plane, 
coloured by, respectively, each point's &e, q P , 9, and /3i value. 
These clearly show that the degeneracy between the slope and Hub- 
ble parameter is tightly connected to the normalisation of the poten- 
tial fee, the axis ratio q p of the isopotential contours, and the source 
position f3, but is completely independent of the orientation angle 
9. Thus, we should not speak just of a r\-h degeneracy, but rather of 
a strong r/-h-bE-q P -/3 degeneracy! This demonstrates that 4-image 
lens systems cannot meaningfully constrain the Hubble parameter 
with significance unless something can be done to at least partially 
break this degeneracy. Constraints on the profile slope or normali- 
sation seem the best hope, as it impossible to learn anything of the 
source position by methods other than lensing. The shape of this 
degeneracy now makes clear the source of the error in Hubble pa- 
rameter and slope returned when the slope is fixed to the wrong 
value, as in Figure [TJ setting the slope is equivalent to drawing a 
line at a particular value of 77 in the scatter plots of Figure [5] and 
taking only the points that fall on that line. If the line is taken at the 
wrong place, completely adequate models will be found that will 
return completely false values for the other model parameters. 

Figure |6]plots the marginalized mean best-fitting and the max- 
imum likelihood values for the Hubble parameter h, the profile 
slope parameter rj, the potential normalisation 6e, and the axis 
ratio of the isopotential ellipses q p obtained for the prolate and 
oblate halos as a function of orientation angle. We find that the 
most-probable values of h and rj vary inversely, with h getting 
smaller and r/ larger with increasing visible ellipticity. The po- 
tential normalisation &e decreases significantly for the prolate ha- 
los and increases significantly for the oblate halos with increas- 
ing ellipticity; this is as expected, because in the maximum ellip- 
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Figure 7. 1-cr (68%) confidence contours for power law model parame- 
ters 77, bE, q P , and /3i with Hubble parameter for a sampling of symmetric 
oblate and prolate halos of triaxial minor axis ratio a = 0.6 at a range of 
rotation angles and thus projected ellipticities. 



ticity orientation prolate halos have minimal mass hidden along 
the line-of-sight, and oblate halos maximum, and the converse in 
the circularly-symmetric orientation at 0°. This result is consistent 
with the behaviour found in previous analyses that weak lensing 
mass estimates are stro ngly effected by unaccounted-for triaxiality 
dCorless & Kind2007l) . The values of the axis ratio of the isopoten- 
tial contours q p exhibit the trend expected, decreasing as the visible 
ellipticity increases from orientation 0° to 90° . 



3.4.2 Hubble parameter as a function oftriaxial orientation: 
posterior probability distributions 

To better understand what is happening to the posterior probabil- 
ity distribution as the triaxial halo orientation is changed, in Figure 
|7] we plot the 1-cr (68%) confidence contours for 77, &g, q p , and 
/3i against h for prolate and oblate halos at several orientations, 
again with visible ellipticity increasing with rotation angle. For 
both oblate and prolate halos, the r\-h contours maintain approx- 
imately the same position, but are tighter at higher visible elliptici- 
ties. The least elliptical case alos shows an extention of its contour 
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to lower values of r\ and higher values of h. The contours for the po- 
tential normalisation bE-h move down to lower normalisations for 
the more elliptical prolate halos, and up to higher normalisations 
for more elliptical oblate halos, as expected given the line-of-sight 
masses of halos of these geometries. Further, they decrease in size 
as ellipticity increases, as with the r\-h contours. For both prolate 
and oblate halos the contours for isopotential axis ratio q p -h move 
down to lower values as the ellipticity increases, as expected. The 
/3i-h contours behave very similarly to those for h-r\, maintaining 
similar positions in parameter space but changing size. 

To understand the positioning and size of these contours we 
consider the nature of the multi-variable degeneracy. It lies in pa- 
rameter space such that high values of the axis ratio q p correspond 
to high values of the slope parameter 77, and low values for the 
Hubble parameter h and source position /3. Now, for halos with 
low amounts of visible ellipticity and thus true values of q p near 
1, there is less space in parameter space for models at values of 
q v higher than the true value, because the posterior distribution is 
truncated with a hard wall of zero probability at q v — 1. Thus, the 
contours will extend proportionally more into the regions of lower 
than true q p , corresponding to lower values of r\ and higher val- 
ues of the Hubble parameter h and source position /3. This trend 
is exactly what is seen as we move from the low-g p halos (rotation 
angles near 90°) to the high-g p halos aligned closer to the circularly 
symmetric (0°) position. 

In addition to their position, the high-g p halos also produce 
larger la contours. This change in contour size cannot be due to 
the lensing strength: in strong lensing, so long as the same num- 
ber of images are resolved, the relative strength of the lens and the 
scale of the image separations does not change the amount of in- 
formation available to constrain the lens system. This is born out 
in the results, in that it is the high-g models in both oblate and pro- 
late cases that exhibit larger contours; in the prolate case these are 
indeed the stronger lenses with higher projected masses, but in the 
oblate case they are the weaker lenses, with very low lensing con- 
vergences. Instead, this again appears to be the natural product of 
the truncation of the posterior probability distribution at q p — 1, 
which shifts and extends the la region towards lower values of 
q p and r\ and to higher values of h. This effect - not an error or 
a problem but a true representation of the shape of the posterior 
probability distribution - is illustrated and explored in more detail 
in Appendix lAl 



3.4.3 Hubble parameter as a function oftriaxial axis ratio 

Again, to understand the importance of the degree of triaxiality, 
now in the general case, we carry out the same analysis as above but 
for prolate and oblate halos of axis ratios a — {0.5, 0.6, 0.7, 0.8}, 
all oriented at 80°. Figure[8]plots the marginalized best-fit parame- 
ters h, r/,bE, and q p as a function of minor triaxial axis ratio a for 
both prolate and oblate halos. Again, increasing the amount of tri- 
axiality does not significantly change the best-fitting slope or Hub- 
ble parameter within the errors. There is a weak trend toward there 
being larger offsets from the true value at smaller visible elliptici- 
ties, the opposite of that seen in the fixed-slope case. This trend can 
be understood as an extension of that seen for this general case as 
a function of orientation angle; again, the la region is shifted to- 
wards lower values of q p and 77 and to higher values of h for lower 
visible ellipticity cases; again see Appendix lAl for further discus- 
sion. 
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Figure 8. Marginalized mean best-fit parameters h, rj, bE, and q p for 
symmetric oblate and prolate lensing halos of minor axis ratio a = 
{0.5, 0.6, 0.7, 0.8} (filled symbols) at a rotation angle of 80°, as well as 
the maximum likelihood parameters (unfilled), plotted as a function of the 
minor axis ratio. The dotted lines show the true parameter values of the 
underlying lensing halo. 
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Figure 9. Marginalized mean best-fitting parameters for symmetric prolate 
halos of minor axis ratio a = 0.6 under three different priors on the density 
profile slope r\, plotted as a function of the rotation angle of the lensing halo: 
at 0° the halos are circularly symmetric in projection, at 90° they show a 
maximum visible ellipticity of 0.6 in the projected isodensity contours. The 
dotted lines show the true parameter values of the underlying lensing halo. 
A tighter prior centred around the true slope value moves all the marginal- 
ized parameter values towards the true value. Such a prior may often be 
justified by independent knowledge of the galaxy density profile. 
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3.4.4 Priors and best-fitting parameters 

These results represent those of the most general analysis, taking 
wide, flat priors on all parameters (0.2 < h < 1.5; 0.5 < r\ < 1.5; 
0.1 < q < 1.0; bE ^ 0; /3 and unconstrained). However, in 
most cases, it is possible and indeed preferable to impose a physi- 
cally realistic prior on the density profile slope, as work on galaxy 
profiles converges to show that nearly all galaxies are very close 
to isothermal. To understand the impact on our results of including 
this physically motivated prior we impose two different Gaussian 
priors on the slope parameter 77, centred on the isothermal value 
■q = 1, with a = 0.1 and a — 0.05. The marginalized best-fit 
parameters found under these prior are shown in Figure [9] As ex- 
pected, the best-fitting parameters migrate towards the true values 
of the system when the prior is tightened on the slope, as the r\ — h 
degeneracy is partially broken in favour of the true values of the 
lensing model employed in the simulations. 

We note a significant offset from the true values in the most- 
probable marginalized parameters under a flat prior; the Hubble pa- 
rameter estimates are uniformly high and the slope estimates uni- 
formly low. Is this a product of triaxiality? Examining the results 
obtained when fitting to halos of different degrees of triaxiality at 
the same orientation (as shown in Figure [8]l and those obtained 
when fits are carried out for fixed slope parameters (as shown in 
Figure QJ, it is clear the answer is no. Most importantly, the values 
of h and r\ obtained in the general model do not change significantly 
with the geometry of the halo - oblate and prolate halos with the 
similar visible axis ratios give very similar results. The visible axis 
ratio has some impact on parameter estimates due to the changing 
shape of the posterior probability distribution with respect to the 
hard wall at q v — 1, but the line-of-sight structure has no impact. 
Thus, the offset in parameters cannot be due to the 3D shape of the 
halo. 

Further, it may seem surprising that the most-probable 
marginalized parameters under a flat prior do not include the true 
values within their la error bars; however, the reliability of the 
result is underlined by the steady movement of the best-fitting pa- 
rameter values to the true underlying values as an ever tighter prior 
is imposed on the profile slope, culminating in the return to the 
true value when a 5-function prior is imposed and the slope set to 
isothermal. Thus, it is not triaxiality, but the large r\-h-bE-q v -fi de- 
generacy, that causes the offset in parameter values. The shapes of 
the parameter contours tell us that, with a flat prior, a much larger 
volume of high-likelihood models reside at high values of h, low 
values of rj, low values of 6b, and low values of q p in parame- 
ter space. The effect is strongest for low projected ellipticity halos 
with values of q v near 1 due to the hard wall at q p = 1, as the con- 
tours are pushed proportionally even more to low values of q v and 
thus lead to even higher best-fitting values of h and lower values 
of r\. Thus the high values of the Hubble parameter are not errors, 
per se, but the true best estimates of the parameter if the degenera- 
cies in the problem are unbroken. It is therefore crucial to break the 
degeneracy in order to obtain meaningful parameter estimates. 



4 CONCLUSIONS 

Both cluster and galaxy halos are predicted to exhibit significant 
triaxiality; in this study we have assessed the impact of neglecting 
this predicted halo shape on Hubble constant estimation in galaxy 
strong lensing analyses, and find it is negligible. Fitting elliptical 
power law models to triaxial oblate and prolate halos using a flex- 
ible MCMC technique shows that, when the density slope of the 



model matches that of the lensing halo, correct Hubble parameter 
values are always recovered within the errors without bias, a result 
that we further explore and explain analytically. In the case of mis- 
matched slopes the presence of projected elliptical structure may 
slightly increase the error in the recovered Hubble parameter values 
as compared to more spherical cases. However, prolate and oblate 
halos of the same projected axis ratio exhibit exactly the same be- 
haviour, in spite of very different line-of-sight geometry, indicating 
that it is increased projected ellipticity and the resultant shifting of 
the posterior towards relatively low values of q v that is significant, 
not 3D triaxial structure. Thus, we find that halo triaxiality can- 
not contribute to the inconsistencies between some lensing-derived 
values of Ho and those of the HST Key Project. Note that we have 
considered very triaxial models motivated by (dark matter) N-body 
simulations, which would represent some of the most triaxial galax- 
ies even when baryonic matter (which pla ys a crucial role on the 
scales relevant to strong lensing, see e.g. JDebattista et al.l2007f) ) is 
included. 

The Hubble parameter as derived by strong lensing is degen- 
erate not only with the slope of the lensing density profile, but also 
with the normalisation of the lensing potential, the projected el- 
lipticity of the lens, and the source position. Accurately constrain- 
ing the profile slope successfully breaks this significant degeneracy. 
However, incorrectly constraining the slope, for example taking an 
average value across a population as true for a given system, will 
lead to well-fitting models with incorrect Hubble parameter values, 
as the degeneracy is broken at the wrong point. This problem is best 
resolved by carrying out statistical studies with many lens systems, 
in which average values of the profile shape are likely to be cor- 
rect when applied across a population. Conversely, if the Hubble 
parameter can be considered well constrained, then strong lensing 
can be successfully used for profile-type determination (i.e. NFW 
vs. Isothermal) without requiring the use of complex triaxial mod- 
els in the fitting process, as our work shows that neither the Hubble 
parameter or density profile slope estimates are biased by neglected 
triaxiality. 

Past studies have concluded that certain effects can exacerbate 
the inconsistencies between lensing-derived Hub ble parame t er val- 
ues and the HST value; for example, Keeton & Zabludoff ( 2004) 
showed that neglect of environment can cause an overestimate in 
Ho values, the opposite trend of that required to account for the low 
values observed in many lensing systems. Thus, the remaining fac- 
tors proposed and still under investigation that may lead to underes- 
timates of the Hubble constant are placed under even greater pres- 
sure to explain the discrepancy. One such promising f acto r is hal o 
substructure, its effects recently demonstrated by Oguril J2007I) ; 
another was highlighted when iDobke. King & FellhaueJ ( 2007) 

showed that, despite findings that in general galaxies strongly tend 

J IT ^ 

towards isothermal profiles ( Koopmans et al. 2006), group interac- 
tions many result in profile steepening that would lead to an under- 
estimate of Ho should an isothermal slope be assumed. Out work 
shows that this effect may be slightly enhanced by highly elliptical 
lenses, but the additional impact will be small. Overall, triaxial- 
ity cannot explain the difficulties encountered in reconciling dis- 
parate lensing re sults for the Hubble parameter with those of other 
techni ques (e.g. lKochanekl r2002. Koc hanekll2003l. l Kochanek et al.1 
2006b). However, the monster degeneracies of the general power 
law model explored in this paper emphasise one reason why achiev- 
ing consistent lensing results is extremely difficult without a good 
external constraint on the density profile slope. 

Large samples of multiply imaged sources with measured time 
delays, deep optical and infrared imaging and spectroscopic data 
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will soon become available. Future space and ground based facili- 
ties such as GAIA, LSST, DES, and further down the line, DUNE 
and the SKA, will greatly enhance our knowledge of the physi- 
cal and environmental properties of lensing galaxies, may generate 
lensing samples large enough to allow cleaner studies using only 
galaxies in isolated environments, and will allow statistical stud- 
ies of the Hubble parameter in large numbers of lensing systems 
(e.g., hundreds of systems in DES) in which the strong degenera- 
cies of the general model parameters can be broken, greatly clari- 
fying galaxy lensing's implications for cosmology. 
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Figure Al. Projections of the posterior probability distribution into the r\-h 
plane of a power law model fit to an isothermal triaxial prolate halo. The 
left-hand panel shows results for lenses of axis ratios (a = 0.6, b = 0.6) 
rotated between orientations of 30° (low ellipticity, q p near 1 ) and 90° (high 
ellipticity, lower q p ); the right-hand panel shows results for lenses of axis 
ratios a = b = {0.5, 0.6, 0.7, 0.8} all oriented at a high-ellipticity orienta- 
tion of 80°, with the low a models having the highest projected ellipticities. 
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Figure A2. Projections of the posterior probability distribution into the r\-h 
plane of a power law model fit to an isothermal triaxial prolate halo with 
three priors on axis ratio q p . The left-hand panel show results for a lens 
with axis ratios (a = 0.6, b = 0.6) at orientation of 90° (high ellipticity, 
lower q p ); the right-hand panel shows results for a lens with axis ratios 
(a = b = 0.5) oriented at a high-ellipticity orientation of 80°. 
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APPENDIX A: PARAMETER ESTIMATES AND 
CONTOUR SHAPES AS A FUNCTION OF PROJECTED 
AXIS RATIO 

When fitting the fully general power-law model, we note a trend 
in which the most-probable parameter values move closer to the 
true values with increasing projected ellipticity. We credit this trend 
to the presence of a hard wall in the prior on the axis ratio q p at 
q p = 1, beyond which the equations governing the lensing poten- 
tial are undefined. In projected halo geometries with low elliptici- 
ties (high q p values), this wall limits the volume of parameter space 
available at values of q p higher than the true value. This is not an 
error; it reflects the true shape of the posterior probability distribu- 
tion. But it does affect the shape and positioning of the resulting 
isoprobability surfaces, pushing them towards lower values of q p 
relative to the true value, and due to the shape of the r;-h-bE-q P -f3 
degeneracy, to higher values of h and /3 and lower values of 77. 

That this effect depends only on the projected shape of the 
potential and not on the 3D geometry is clear when the contours 
of prolate halos of one geometry projected at different orientation 
angles are compared to the contours of prolate halos of different 
underlying geometries (a range of minor axis ratios) all projected 
at the same orientation. The r) — h plane la contours for these two 
cases are shown in Figure lATl In both cases, the projected ellipticity 
varies from high to low, but the underlying geometry leading to that 
change is fundamentally different between the two. However, the 
exact same trend is visible, with the low ellipticity, high-%, cases 
exhibiting larger contours extended toward high values of h and 
low values of rj. 

To further confirm that this behaviour is indeed attributable 
to the necessary imposition of a hard wall of zero probability at 
q p = 1, we look for the opposite effect in the high ellipticity, low 
q p cases, by artificially imposing a hard wall of zero probability 
at q p — {0.8, 0.82}. If our explanation is correct, the impositions 
of these priors should produce a similar but opposite effect in the 
resulting 7] — h contours, shifting them to lower values of h and 
higher values of r\. The results for both the constant geometry and 
constant orientation cases are shown in Figure fA2] and indeed the 
contours behave as expected, with the contours moving more and 
more to higher values of r\ as the prior is tightened. 

Thus, the trends seen when lensing through both halos of the 
same geometries oriented to exhibit different projected ellipticities 
(Figures[6]&[7J and of different geometries oriented similarly (Fig- 
ure [TJ - that systems with projected axis ratios close to q p = 1 re- 
turn parameter values further away from the true underlying value 
of the system than do their more elliptical counterparts - is at- 
tributable to the real shape of the posterior probability distribution 
and is a function only of projected ellipticity and independent of 
3D structure. 



